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ABSTRACT 

I develop a model for the abundance and spatial distribution of dark matter subclumps. 
The model shows that subclumps of massive parent halos formed at earlier times than 
subclumps of the same mass in lower mass parents; equivalently, halos in dense regions 
at a given time formed earlier than halos of the same mass in less dense regions. This 
may provide the basis for interpreting recent observations which indicate that the 
stellar populations of the most massive elliptical galaxies are also the oldest. 
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1 INTRODUCTION 

In hierarchical clustering models about ten percent of the 
mass of a dark matter halo is expected to be in the form 
of subclumps which are the remnants of objects which fell 
in and were tidally stripped as the halo was assembled. 
Whereas models of the abundance of halos as a function 
of mass are quite well developed (Press & Schechter 1974; 
Sheth & Tormen 1999; Jenkins et al. 2001; Sheth, Mo & Tor- 
men 2001), models of the subclump abundances are not. The 
main goal of this paper is to derive a model for the distri- 
bution of subclump masses. Since it is likely that such sub- 
clumps are associated with sites of galaxy formation (White 
& Rees 1978; Colin et al. 1999) such a model is a useful 
step towards interpretting galaxy luminosity functions. The 
model is developed in Sections 2 and 3, and Section 4 illus- 
trates the results. The model presented here may provide 
a useful framework within which to discuss the luminosity 
function of galaxies, as well as the formation of massive ellip- 
tical galaxies. These possibilities are discussed in Section 5, 
which also highlights shortcomings of the model, and out- 
lines possible improvements. 



2 THE POISSON MODEL 

Sheth & Pitman (1997) described a model in which the num- 
ber density of clumps containing m particles is 



n{m,b) = (1 - b) 



{mby 



(1) 



where b which grows from an initial value of zero towards 
an upper limit of unity, plays the role of a time variable. 
They showed that the rate of change of the number density 



of m-halos could be thought of as the difference between 
formation and destruction rates: 
dn{m,b) _ dnform(6!m-) dndcst(b\m) 



db 
where 

dnfoi.m(b|m) 
db 

dndest(b|m) 



db 



db 



(m - 



■ n{m, b) 



n{m, b). 



(2) 

(3) 
(4) 



db V 1-b/ 

In this model, an object of mass M at some late time bo was 
previously (i.e., at some epoch bi < bo) made up of many 
smaller pieces. The number of m-subclumps of an i\f-parent 
is 

where B = bi/bo < 1. Extending their model to compute 
the formation and destruction rates of these m-subclumps 
yields: 

(m-1) 



(5) 



dNf^rmiB\m, M) 
dB 

djVdo.t(g|m,M) 
dB 



B 

(1 ^ m/M) 
(l-mB/M) 



N{m\M, B) 



(6) 
(7) 

N{'m\M,B). (8) 



1 - B 

Details are in the Appendix. In the limit that m <^ M, equa- 
tion (8) reduces to the unconstrained rates. As a consistency 
check, note that 



diVform(-B|m, Af) dnform(bi|m) 
^ -n(m, bo) = 



^ bodS 

M — m 



dbi 



(9) 
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where we have used the fact that dB — dbi/bo. This shows 
that the number of m-subclumps which form within Ap- 
parent halos during the interval dbi around 61, times the 
density of M-parcnts, when summed over all M does indeed 
give the total formation rate density of m-clumps. A similar 
expression holds for the destruction rates. 

In this model, the total number density of m halos which 
will ever form is 



n(m) 



db 





(m 



dnform(&|?w) 

d6 



r(m - 1,0, m) 



r(m, 0,7 



(10) 



where the expressions in square brackets denote incomplete 
gamma functions, with parameters m — 1 and m respec- 
tively. In the limit of large m, this total number reduces to 
1/V27rm3. 

Similarly, integrating dNform/dB over B gives the to- 
tal number of m-subclumps which ever formed within an 
M-halo. With suitable choices of m and M, the value of 
this integral can be used to estimate the typical number 
of galaxy-sized subclumps which are today within virialized 
clusters. Notice that this total number is independent of 60; 
i.e., it depends only on the mass M, and not on the time at 
which the parent M-halo formed. 



3 THE CONTINUUM LIMIT 

Setting b = 1/(1 -I- (5ac), and taking m ^ 1 and 5ac <C 1 
makes n(m, b) resemble the Press-Schechter (1974) model of 
the halo mass function. Halo abundances in this model are 
usually written as 



mn{m\t) , dcr (m) 6sc{t) aim) 
_ dm = — 57 — r ^= exp 

p cr^(m) ^/2^ 



Slit) 



(11) 



where p is the comoving mass density, 5sc(^) is the linear 
theory overdensity required for spherical collapse at t (it is 
1.686 {to/t)'^^^ in an Einstein-de Sitter model). 



cr^(m) 



= 4,/ 



Ah, 

^k'P{k)\W{kR)f 



(12) 



is the linear theory variance of the density field at the present 
time to when smoothed with a tophat filter on the scale 
R = (3m/47^p)l/^ and P{k) oc k", with n = for the 
Poisson model. The use of the tophat filter makes W{x) = 
(3/a;'^)[sin(a;) — a; cos(a::)]. 

The usual Press-Schechter argument is to assume that 
expressions derived from Poisson or white-noise initial con- 
ditions are also valid for arbitrary Gaussian fluctuation flelds 
(i.e., P{k) may depend on k). Therefore, where possible in 
what follows, we will write all expressions in variables which 
are independent of the power spectrum. 

The continuum limit of equation (5) shows that 
mN{m, t\M, T) dm/M is given by an expression like the one 
above, but with the replacements cr^irn) — a'^{M) and 

Sac — > Ssc(t) — Ssc{T). This is the same as the conditional 
mass function of Bower (1990), Bond et al. (1991) and Lacey 
& Cole (1993). The continuum limits of the formation time 
distributions given above are 



dnform(t|m) 
dt 



dt — > mn{m\t) 



dt 



(13) 



and 

dNio,^(t\m,M,T) 
dt 



dt mN{m,t\M, T) 



dt, 



(14) 



where m < M and t <T. Furthermore, one can verify that 

(15) 



dNforrn{t\m,M,T) dnforin(t|m) 
— - n{m, T) dm = 



dt 



dt 



Expressions similar to equations (13) and (14) above for the 
formation time distribution have been derived by Percival & 
Miller (1999) from a quite difiierent approach. However, they 
were only able to describe the shape of the formation time 
distribution, not the normalization. What we have shown 
here is that the normalization can be derived from the earlier 
results of Sheth & Pitman (1997). In hindsight, notice that 
requiring the identity above (equation 15) to be satisfied also 
provides a way to determine the normalization factors. 

Note that these expressions for the formation times 
show that subclumps of the same mass form at earlier times 
in massive halos than in less massive ones. We will return 
to this shortly. 



3.1 Subclump abundances 

The total number of m subclumps which form within an 
M-halo is 



7V,ub(m|M,T) 



Jo 



dt 



dNtorm{t\m,M, T) 

dt 



(16) 



If these subclumps survive the processes of tidal stripping 

and dynamical friction as they merge into their parent halo 
(an unrealistic approximation!), then the total number den- 
sity of m-subclumps at time T is 



/•oo 

,(m,r)= / dMn(M,T)iVaub(m|M,T). (17) 



^sub 



3.2 Large-scale spatial distribution of subclumps 

The distribution of M-halos is biased with respect to the 
overall distribution of dark matter. On scales which are large 
compared to the typical diameter of a halo, this bias is fac- 
tor is approximately independent of scale, and is usually 
denoted b{M,T). If there are Nsuh{m\M,T) subclumps for 
each M-halo, then, on large scales, the bias factor associated 
with the subclump distribution is given by a simple counting 
argument: 

, , _ i::dMn{M)NsMm\M,T)b{M) 
6.„b(m,i)- joo^^^^^^r^^^^^^^^^^^r^^ ■ (1«) 

The subclump distribution on smaller scales can be esti- 
mated if one assumes a model for the density run of sub- 
clumps around their parent halos. For instance, one might 
assume that they trace the density profile of the mass, or 
that this distribution is a non-singular isothermal sphere. 
Once this density profile has been specified, the small scale 
distribution of the subclumps can be described rather com- 
pactly using the halo model of large scale structure (see the 
recent review in Cooray & Sheth 2002). 
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3.3 Dependence on local density 

It is sometimes of interest to estimate the distribution of 
subclumps as a function of local density. We develop two 
approximations for this quantity. The first is a simple count- 
ing argument which exploits the fact that the number of m- 
subclumps of an Af-parcnt docs not depend on any quanti- 
ties other than m and M. If the local ovcrdcnsity 5 is defined 
on sufficiently large scales that Ms = pV{l + S) is much 
larger than the mass of any parent halo, then the density of 
subclumps in such regions is 



nsub{m,T\6) = 



Ms 



dM n(M, T\5) Af,ub(m|M, T) 



(19) 



where the density dependent parent halo mass function 
n{M,T\5) can be estimated following Mo & White (1996), 
Kauffmann & Lemson (1999) and Sheth & Tormen (2002). 
On very large scales. Ma 3> M , the parent halo mass func- 
tion reduces to n{M,T\5) -> [I + h{M)5] n{M,T). Thus, on 
these scales, nsnb{m,T\S) [1 + bsub{m,T)5]nsuh{M,T), 
with 6sub(m, T) given by the expression derived in the pre- 
vious subsection. On intermediate scales, this simple ex- 
pression for n{M,T\6) is no longer valid, but equation (19) 
will remain accurate provided that the full expression for 
n{M,T\S) is used. 

Equation (19) becomes inaccurate once the cell size V 
on which 6 is evaluated is comparable to the diameter of a 
typical halo. On such small scales, some cells may only con- 
tain portions of parent halos, and hence the simple count- 
ing argument is no longer accurate. Then it is more accu- 
rate to simply compute NsubifnlMgjTs), where Ts is ob- 
tained by using the spherical collapse relation for the linear 
and nonlinear ovcrdonsitics, 5o{5) (Mo & White 1996 give 
a convenient fitting formula for this relation), and then set 
Ssc{Ts) = 5o{d). Note that, in contrast to equation (16), here 
the upper limit of the integral, T, is not the same as the time 
when the parent halo is identified, Tg. This is because not 
all the subclumps which will form within Mg by Tg have 
formed at T < T^. Thus, 



N,^b{m\Ms,Ts) 



f 

Jo 



dt 



dNfo,^(t\m,Ms,Ts) 
dt 



(20) 



Dividing this quantity by nsub{'m,T)V provides an estimate 
of how biased the subclump distribution in dense cells is rel- 
ative to the average. Writing the density- dependent abun- 
dance this way shows explicitly that subhalos of a given mass 
in dense regions formed at earlier times than did subhalos 
of the same mass in less dense regions. 



4 EXPLICIT RESULTS 

For the parent halo mass function in equation (11), all the 
integrals above can be performed analytically: 



N^ub{m\M,T) dm 



M 



da^{m) 



2-K ^a'^m)-a'^{M) 



(21) 



(note that this is independent of T) and so the abundance 
of subclumps is 

sUt) 



nsub{m,T) , dfT^(m) 
-dm= — exp 
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(22) 




Figure 1 . Large-scale bias factor of the parent halos (dashed line, 
from Sheth & Tormen 1999) and for subclumps (solid line). The 
lower and upper dotted curves show the parent and subclump 
bias factors associated with equations (23) and (24). 



The halo bias factor associated with equation (11) is 



b{M,T) = l + 



5l(T)/a\M) 



(23) 



(Cole & Kaiser 1989; Mo & White 1996; Sheth & Tormen 

1999), so the large scale bias factor of the subclumps, com- 
puted using the simple counting argument of equation (18) 
is 



6,ub(m,T) = 1-1- 



(24) 



Whereas the halo bias factor can be less than unity, the bias 
factor of the subclumps is always greater than unity. The 
more accurate model for subclump biasing (equation 20), 
which should be accurate on smaller scales, can also be 
solved analytically: 



N,ubim\Ms,Ts) 
n,ub{m,T)V 



(1 + 5) 



V'l-o-2(Mi)/cr2(m) 



exp[-i.|/2) 



where vg = 



exp[-5sc(T)V2a2(m)]' 

2 



(25) 



5sc{T)-5o{5) 



0-2 (m) - (T^iMs) 



If this ratio is expressed as a Taylor series in S, then the 
coefficient of the term which is proportional to d is the same 
as the linear bias factor derived from the counting argument 
(equation 24). To sec this explicitly, consider the limit of 
large cells, for which Ms is large simply because V is large. 
Then a{Ms)/(7{m) — » and 5o{S) ~ 5 < (5sc, so that 



N,^i,im\Ms,Ts) 
n^ub{m,T)V 



(l + <5) 



l + <5 



(26) 
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Numerical simulations show that the spherical collapse 
based equation (11) is a good but not perfect description of 
the number density of parent halos. A more accurate formula 
is 

rnn{rn\t) 



dm ■ 




-^U[i + (« 



-2p 



](27) 



(Sheth & Tormen 1999), where v = 4c(i)/o-(m), a « 0.71, 
p = 0.3 and A = 0.322 insures that the distribution is nor- 
malized to unity. Sheth, Mo & Tormen (2001) argue that 
this expression may be related to models in which halos 
form from an ellipsoidal collapse. Percival, Miller & Pea- 
cock (2000) show that the insertion of equation (27) in equa- 
tion (13) provides a better description of halo formation in 
simulations than does equation (11). Therefore, it probably 
provides a more accurate model of the subclump abundances 
and bias factors as well. To use this requires a model for the 
analog of the conditional mass function (equation 5); Sheth 
& Tormen (2002) argue that simply changing variables in 
n{m,t), as is appropriate when equation (11) is the mass 
function, while incorrect, is not a bad approximation. In this 
approximation, the analog of equation (21) is again propor- 
tional to M I ^ a'^(m) — CT^(M), the subclump mass function 
(integrated over parent masses M) can be written in terms 
of hypergeometric functions, as can the large-scale bias fac- 
tor. The expressions arc lengthy, so I haven't reproduced 
them here. Figure 1 compares the large-scale bias factor of 
the parent halos (dashed, from Sheth & Tormen 1999), with 
that for the subclumps (solid). The lower and upper dotted 
curves show the parent and subclump bias factors associated 
with equations (23) and (24). 
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Figure 2. Distribution of formation times for subclumps of mass 
m/M* = 0.01 for various choices of the parent mass Af at z = 0. A 
power-law spectrum with n = —2 was used; lookback times were 
computed assuming a flat cosmology with CIq = 0.3, Aq = 1 — 
and Hubble constant Hq = 70 km Mpc"'^. To set the scale, 
the value of M« in such a ACDM model would be ~ 10^^h~^MQ. 



5 DISCUSSION 

I derived a model of the subclump distribution under the as- 
sumption that there are no processes by which a subclump 
can lose mass as it falls in to a larger system. Although this 
idealization is not realistic, it should be thought of as provid- 
ing the initial conditions for more sophisticated calculations 
which do incorporate the effects of tidal stripping and dy- 
namical friction. For instance, the cumulative distribution 
of subclumps in this model is 

f'^ / 2Ma{m) f a^{M) 
/ Afeub m M, T dm = W 1 - ( . 

When m -C M, this scales as Ma{m), which is rather dif- 
ferent from the scaling seen in simulations after tidal 
stripping. Nevertheless, the model is useful because it pro- 
vides simple closed-form expressions for how the abundances 
and spatial distributions of the subclumps differ from those 
of their parent halos. 

If mergers are unimportant, then equation (22) can be 
turned into an estimate of the luminosity function. This 
can be done either by assuming a mean mass-to-light ra- 
tio, or by assuming some relation for how the light-curve 
evolves with time, and then convolving with the formation 
time distributions derived here. I have not pursued this fur- 
ther because this model for the subclump distribution is 
not sufficiently realistic. For example, J dmmnsnb{m)/ p is 
generally greater than unity. This is a consequence of the 
fact that objects counted as having formed with mass m 



may previously have been counted as having formed with 
mass m' < m. A more realistic model of the subclump dis- 
tribution would avoid this double counting: one attractive 
model in this regard is to study how the mass of the most 
massive progenitor subclump decreases with lookback time. 
One then labels as subclumps all the objects which merge 
with the most massive progenitor. This distribution of sub- 
clumps can bo described following results in Nusser & Sheth 
(1999) and is the subject of work in progress. 

There is one respect in which the model developed here 
is realistic: it shows that subclumps of massive parent halos 
formed at earlier times than subclumps of the same mass in 
lower mass parents, or equivalently, that halos in dense re- 
gions at a given time formed earlier than halos of the same 
mass in less dense regions. This has an interesting conse- 
quence for the following simple model of elliptical galaxies. 
Suppose that gas can only cool and form stars within halos 
which are more massive than some minimum mass m. If this 
happened at high redshift, then the halos within which the 
stars formed will have subsequently merged with other ha- 
los to make more massive objects. If wo treat the separate 
parcels of stars as representing the subclumps m, then we 
have a model in which massive galaxies contain older stars. 
Moreover, since higher redshifts correspond to smaller in- 
tervals of time, the distribution of formation times will be 
narrower for the stars which form in more massive halos, 
and broader for the stars which form in lower mass halos; 
Figure 2 illustrates. This is in qualitative agreement with 
recent work (Thomas, Maraston & Bender 2002). 
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term is the number of subclumps containing rn — j particles, 
given that there are only M j remaining particles from 
which they could have been formed (see Sheth & Lemson 
1999 for a detailed discussion); this is why the second term 
is not N{m—j, bi\M, bo). For similar reasons, the destruction 
rate is obtained from 

M — m—l 

N{m,bi\M,bo) {m+j)N{j,bi\M-m,b') 

3 = 1 

rather than N{m, bi\M, bo) J2^J[""~\m + j)N{j, bi\M, bo). 



APPENDIX A: THE ADDITIVE 
COAGULATION MODEL 

In the model discussed in the main text, halos grow by bi- 
nary mergers with one another. The formation rate of a halo 
with M particles is given by evaluating 

— n(m, b)n{M — m,b), 

m=l 

which expresses the fact that formation in this model hap- 
pens by binary mergers, with probability of merger propor- 
tional to the mass of the merging pair, and the abundance 
of each type of clump. The associated destruction rate is 

related to 

oo 

n{M,b)y^^{M + m)n{m,b). 

m=l 

This system of equations for n{M, b) is solved by equa- 
tion (1); this was one of the results in Sheth & Pitman 
(1997). 

The conditional formation rate of m-subclumps within 
an M-parent, while given by a similar expression, is slightly 
different. It is 

M — m— 1 

^ N{j,bi\M,bo)N{m-j,bi\M-j,b'), 

3 = 1 

where 1/6' = (M - j)/{Mbo - jbi). The first term is the 
number of j-subclumps within the M-halo, and the second 
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